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ABSTRACT 

(Ultra) Luminous Infrared Galaxies ((U)LIRGs) are objects characterized by their extreme infrared 
(8-1000 p m) luminosities ( Lurg > 10 n L o and Luurg > 10 12 L©). The Herschel Comprehensive 
ULIRG Emission Survey (HerCULES; PI van der Werf) presents a representative flux-limited sample 
of 29 (U)LIRGs that spans the full luminosity range of these objects (10 11 < L© > 10 13 ). With the Her¬ 
schel Space Observatory, we observe [CII] 157 pm, [01] 63 pm, and [OI] 145 pm line emission with 
PACS, CO J=4-3 through J=13-12, [Cl] 370 pm, and [Cl] 609 pm with SPIRE, and low-J CO transitions 
with ground-based telescopes. The CO ladders of the sample are separated into three classes based on 
their excitation level. In 13 of the galaxies, the [OI] 63 pm emission line is self absorbed. Comparing the 
CO excitation to the IRAS 60/100 pm ratio and to far infrared luminosity, we find that the CO excitation 
is more correlated to the far infrared colors. We present cooling budgets for the galaxies and find fine- 
structure line flux deficits in the [CII], [Sill], [OI], and [Cl] lines in the objects with the highest far IR 
fluxes, but do not observe this for CO 4 < J upp < 13. In order to study the heating of the molecular gas, 
we present a combination of three diagnostic quantities to help determine the dominant heating source. 
Using the CO excitation, the CO J=l-0 linewidth, and the AGN contribution, we conclude that galax¬ 
ies with large CO linewidths always have high-excitation CO ladders, and often low AGN contributions, 
suggesting that mechanical heating is important. 

Subject headings: 


Leiden Observatory, Leiden University, P.O. Box 9513, 2300 
RA Leiden, The Netherlands 

department of Earth and Space Sciences, Chalmers University 
of Technology, Onsala Observatory, 43994 Onsala, Sweden 

3 Spitzer Science Center, California Institute of Technology, MS 
220-6, Pasadena, CA 91125, USA 

institute for Astronomy, Astrophysics, Space Applications & 
Remote Sensing, National Observatory of Athens, P Penteli, 15236 
Athens, Greece 

5 Department of Astronomy, University of Virginia, PO. Box 
400325, Charlottesville, VA 22904, USA 

6 Naval Research Laboratory, Remote Sensing Division, 4555 
Overlook Ave SW, Washington, DC 20375, USA 

7 Purple Mountain Observatory, Chinese Academy of Sciences 
(CAS), 2 West Beijing Road, Nanjing 210008, China 

8 Universidad de Alcala, Departamento de Fsica y Matematicas, 
Campus Universitario, 28871 Alcala de Henares, Madrid, Spain 
department of Physics and Astronomy, University College Lon¬ 


don, Gower Street, London WC1E 6BT, UK 

^Department of Astronomy, University of Maryland, College 
Park, MD 20742, USA 

11 Max-Planck-Institut fur Radioastronomie, Auf dem Hiigel 16, 
Bonn, D-53121, Germany 

12 Scientific Support Office, ESTEC/SRE-S,Keplerlaan 1, NL- 
2201 AZ, Noordwijk, The Netherlands 

13 Instituto Radioastronorma Milimetrica (IRAM), Av. Divina 
Pastora 7, Nucleo Central, 18012 Granada, Spain 

14 Institute for Space Imaging Science, Department of Physics and 
Astronomy, University of Lethbridge, Lethbridge, AB, Canada 
15 Institute for Astronomy, 2680 Woodlawn Drive, University of 
Hawaii, Honolulu, HI 96822 

16 Harvard-Smithsonian Center for Astrohpysics, 60 Garden 
Street, Cambridge, MA 02138, USA 
17 Kapteyn Astronomical Institute, P.O. Box 800, NL-9700 AV 
Groningen, The Netherlands 


1 



1. Introduction 

(Ultra) Luminous Infrared Galaxies ((U)LIRGs) in 
the local universe are remarkable galaxies exhibiting 
an extremely high infrared luminosity, L« looopim > 
10 n L o for LIRGs and > 10 12 L G for 

ULIRGs. Luminous infrared galaxies were first iden¬ 
tified in large numbers with observations from the 
InfraRed Astronomical Satellite (IRAS), which was 
launched in 1983 ( |Houck et al.||1985| ). After the dis¬ 
covery that these objects all contain massive amounts 
of molecular gas (Sanders et al. |1988a Veilleux et al. 
2002), detailed studies of the spectroscopic cooling 
lines were carried out with the Infrared Space Ob¬ 
servatory (ISO;|Malhotra et al.|(jl997ji; Luhman et al. 
( fl998l >; |HeIbu et al.| ( |2001| >; |Malhotra et al.| ( |2001 [ >T 
Luhman et al. ( 2003| l; Abel et ah] (|2009|l), ground 


1996[ Veilleux et al. 2002| l. Regardless of the 


based observations of [Cl] (Gerin & Phillips 


2000), Spitzer Space Telescope (Armus et al. 


1998 


2009 


Dfaz-Santos et al 7| |20 11 [ jStierwalt et al.||2013| ) and 


the Herschel Space Observatory (|Gracia-Carpio et al. 


|201 1[ |Dfaz-Santos et al.|20131 jFarrah et al.H201 3] ' Lu 

|et al.|2014t|Dfaz-Santos et al.|2014 ). In the local uni¬ 
verse ULIRGs are rare (Soifer & Neugebauer 1991] ), 
but at higher redshifts (z > 1) they represent most of 
the cosmic infrared background and are the dominant 
source of star formation up to z=2 (jCaputi et al.|2007 


Magnelli et al. 2011^ Berta et al. |2011a|b[ Magnelli 


et al. 2013[ Gruppioni et al. 20 13| )~ Locally, these 


objects are hosts to intense starbursts, and/or active 
galactic nuclei (AGN), and often are part of a merg 
ing galaxy group (Armus et al. 1987 jSanders et al. 


1988b| Barnes & Hemquist|1992 Sanders & Mirabel 
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ous heating processes available, however, the luminos¬ 
ity of most local (U)LIRGs seem to be energetically 
driven by starbursts (Genzel et al. 1998[ Downes & 


|Solomon|[l998l IVeilleux et al.||1999| |2002[ |Gao & 


|Solomon|[2004l |Veilleux et ak1|2009| ). (U)LIRGs are 
also thought to represent the transitional phase in evo¬ 
lution from a starburst galaxy to elliptical/lenticular 


galaxies (Sanders et al. 

1988b Genzel et al. 2001 

Tac- 

coni et al. 

2002 Rothberg & Fischer 2010, Rot 

hberg 

et al. 2013 

), and thus must quench their star formation 


during this period. In fact, some evidence for this was 
found in the discovery of massive molecular outflows 


with the Herschel Space Observatory (Fischer et al. 


20|0] Sturm et a l.||2011 [ Spoon et al. j2013 Veilleux 


et a l.[2013 Gonzalez-Alfonso et al.|2014| as well as as 

by ground-based telescopes (e.g., Feruglio et a l.|2010 
|Weifi et al.|2dl2] ). 

Since (U)LIRGs offer a unique insight into this 
transitional phase from star-forming to quiescent 
galaxies, understanding which mechanisms are affect¬ 
ing the star-forming gas is crucial. Many studies of 
the star forming gas in (U)LIRGs have been made 
since its universal presence in (U)LIRGs was deter- 


mined (Sanders et al. 

1991 1 Sanders & Mirabel 1996; 

Solomon et al. 1997 

). In general, gas is heated by 


either radiation (i.e. UV photons. X-ray photons), 
energetic particles (cosmic rays) or mechanical pro¬ 
cesses (i.e. turbulence, stellar winds, outflows, super¬ 
novae). The interplay between these heating sources 
can account for the extreme environments found in 
(U)LIRGs, in comparison to less intense star form¬ 
ing environments ( | Aalto et al.|1991 1995a). The high 
amount of energy injected into the gas in these galaxies 
is displayed by emission lines that serve as a coolant 
along with infrared dust emission. The emission lines 
responsible for most of the gas cooling are the [CII] 
line at 157 pm ( 2 Pi /2 - 1 P\/i), the [OI] line at 63 pm 
( 3 P i - 3 Pi), and CO (rotational transitions). The Her¬ 
schel Space Observatory has, for the first time, pro¬ 
vided astronomers with simultaneous access to these 
important far infrared cooling lines and the CO rota¬ 
tional ladder (CO ladder) in (U)LIRGs. Using the mul¬ 
tiple rotational transitions of CO from J=l-0 through 
J=13-12, the density, temperature, column density, 
and mass (with the addition of 13 CO) can be esti¬ 


mated (eg. Rangwala et al. 2011 Spinoglio et al. 2012 


Rigo poulou et al.|2013| Papadopoulos et al.|2014| ). In 

some cases, it is possible to even discern specifically 
the heating mechanism (|Loenen et al.||2008] |Hailey- 
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Dunsheath et al.|2008[ van der Werf et al. 2010a||Mei-| 

jerink et al. 2013; Lu et al. 2014 

Rosenberg et al. 

2014a Pereira-Santaella et al. 2014 
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In this paper, we introduce observations of all major 
neutral gas cooling line of a representative sample, the 
HerCULES sample, of local (U)LIRGs spanning the 
luminosity range from 10 11 <L fir < 10 13 L G . In Sec¬ 
tion [2j we present the HerCULES sample and obser¬ 
vations from the Herschel/SPIKE and HerschelfPACS 
spectrometers, which include [CII], [OI] 63/zm, [OI] 
145/rm, CO (4 < J upp < 13), and [Cl] 370 pm and 609 
pm, using the cosmological parameters Ho=70 km s ' 1 
Mpc, Q. vaccum =0.12, and Q. matter =0.28. In this paper we 
focus on the main neutral gas cooling lines. We there¬ 
fore do not analyze the [Nil] lines, which arise in ion¬ 
ized gas, or other molecular lines which do not affect 
the thermal balance. Specifically, we do not discuss 
H 2 O since in the cases where these lines are bright, 
there is strong evidence that they are radiatively ex¬ 
cited ([Gonzalez-Alfonso et al.|20l0||van der Werf et al.| 


2011 11 Yang et al.|2013 ! and do not remove kinetic en¬ 
ergy from the gas and thus do not contribute to the 
cooling. We show spectra for three sample galaxies 
that represent three different classes of excitation, and 
the CO ladders for the full sample in Section [3] Using 
the full sample, in Section [4] we analyze the gas ex¬ 
citation, cooling budget of the sample, and a diagram 
for determining additional heating mechanisms for the 
gas. Our conclusions are presented in Section [5] 


2. Observations 


a range of objects including starburst galaxies, AGNs, 
and composite sources, and covering also a range of 
IRAS 60/100 micron ratios. The full list of included 
galaxies and their respective properties can be found 
in Table [T] The infrared luminosity and the luminosity 
distance are from |Armus et al.| ( [2009j >. 

In order to obtain a comprehensive view of the CO 
emission and the cooling budget of these galaxies, 
we proposed Herschel/ SPIRE spectroscopy (for the 
CO ladder) and HerschelfPACS spectroscopy (for the 
[CII] and [OI] fine structure lines) of the entire sam¬ 
ple, unless PACS observations were already observed 
as part of another program. In addition to the galax¬ 
ies observed for HerCULES, we have included NGC 
4418, NGC 1068, and Arp 220 for completeness. This 
project was approved as a Key Project on the Herschel 
Space Observatory, under the name Herschel Compre¬ 
hensive (U)LIRG Emission Survey (HerCULES - PI. 
Van der Werf). Key elements of HerCULES are: 

• a representative flux-limited sample of local 
LIRGs and ULIRGs; 

• comprehensive coverage of the SPIRE spectral 
range at the highest spectral resolution mode 
(covering the CO ladder, [Cl] and [Nil] fine 
structure lines, and any other bright features 
such as H 2 O lines); 

• comprehensive coverage of the key fine-structure 
cooling lines [CII] and [OI] with PACS observa¬ 
tions; 


2.1. The HerCULES sample 

The sample was chosen from the IRAS Revised 
Bright Galaxy Sample (RBGS), which contains all 629 
extragalactic sources with IRAS 60 pm flux density 
S 60 > 5.24 Jy in the (IRAS-covered) sky at Galac¬ 
tic latitudes \h\ > 5 ( |Sanders et al.|[2003l ). From the 
IRAS RBGS we select a sub-sample applying limits 
both in Sgo and L :r : at luminosities : L IR > 10 12 L 0 
(ULIRGs), all sources with S 60 > 11.65 Jy are in¬ 
cluded, while at luminosities 10 11 L Q < Li R < 10 12 
L 0 (LIRGs), sources with Sgo > 16.4 Jy are in¬ 
cluded. From this flux-limited representative parent 
sample of 32 targets, we removed three LIRGs for 
which no ground-based CO data are available, with the 
exception of ESO 173-G015, IRAS 13120-5453, and 
MCG+12-02-001. The resulting representative flux- 
limited sample consists of 21 LIRGs and 8 ULIRGs. 
The sample covers a factor of 32 in L IR and contains 


Details about the galaxy type and observation ID 
can be seen in Table Q] We have included ob¬ 
servations from other programs (KPGT_esturm_lK, 
KPGT_cwilso01_l, OTl_larmus_l, OTl_shaileyd_l) 
to help realize the complete flux-limited sample. The 
references for these observations are also in Table |T] 


2.2. Herschel/SPIRE observations 


Spectra were obtained with the Spectral and Photo¬ 
metric Imaging Receiver and Fourier-Transform Spec¬ 
trometer (SPIRE-FTS, |Griffin et al.|2010 1 on board the 
Herschel Space Observatory (Pilbratt et al. |2010| for 
the full HerCULES sample. The observations were 
carried out in staring mode with the galaxy nucleus on 
the central pixel of the detector array, with a beam size 
varying from 17”-42” for the CO transitions. The high 
spectral resolution mode was used with a resolution of 
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1.2 GHz over the two observing bands. The low fre¬ 
quency focal plane array (Long Wavelength Spectrom¬ 
eter Array, SLW) covers v=447-989 GHz (4=671-303 
H m) and the high frequency focal plane array (Short 
Wavelength Spectrometer Array, SSW) covers v=958- 
1545 GHz (4=313-194 jum), and together they include 
the CO / = 4 - 3 to CO J - 13-12 lines. For galaxies 
with z > 0.03, the rest frequency of the J=4-3 tran¬ 
sition falls short of the SPIRE coverage. All galaxies 
were observed in the sparse observing mode besides 
NGC 4418, which was observed in the intermediate 
mode. 


The data were reduced using version 13.0 of the 
Herschel Interactive Processing Environment (HIPE). 
Initial processing steps included timeline deglitching, 
linearity correction, clipping of saturated points, time- 
domain phase correction, and inteferogram baseline 
subtraction. After a second deglitching step and in- 
terferogram phase correction, the interferograms were 
Fourier transformed, and the thermal emission from 
instrument and telescope was removed from the result¬ 
ing spectra. The averaged spectra were flux calibrated 
as point sources using the calibration tree associated 
with HIPE 13.0. Following these steps a ’’dark” spec¬ 
trum was subtracted, to remove any residual emission 
from the telescope and the instrument. Since the emis¬ 
sion of most of our sources is contained entirely in the 
central pixel of the detector arrays, a ’’dark” spectrum 
was constructed by spectrally smoothing a combina¬ 
tion of several off-axis pixels. For extended targets, 
where the off-axis pixels contain emission, the dark 
was obtained from a deep blank-sky observation ob¬ 
tained on the same observing day. We compared the 
two methods and found no differences, but the noise 
was smaller using the smooth off axis pixels, in the 
case of the compact targets. 


For all extended sources (Arp 299, ESO 173- 
G015,MCG+12-02-001, Mrk 331, NGC 1068, NGC 
1365, NGC 2146, NGC 3256, NGC 5135, and NGC 


7771), an aperture correction is necessary to compen¬ 


sate for the wavelength dependent beam size (Makiwa 


et al. 2013|). We defined a source as extended using 


LABOCA or SCUBA 350 or 450 gm (respectively) 
maps with 8” resolution. We convolved the 8” resolu¬ 
tion maps with the SPIRE FTS resolution, and if the 
galaxy was more extended than the smallest SPIRE 
beam size, we defined it as extended. In order to 
correct for the extended nature of these sources, we 
employ HIPE’s semiExtendedCorrector tool (SECT). 
This tool ’derives’ an intrinsic source size by iterating 


over different source sizes until it finds one that pro¬ 
vides a good match in the overlap range of the two ob¬ 
serving bands near 1000 GHz, and is further discussed 
in [Wu et al.| ( 12013} . We set the Gaussian reference 
beam to 42”, the largest SPIRE beamsize. The beam- 
size corrected flux values for the 10 extended sources 
are listed in Table [2] along with the compact sources. 
We note that the error in the extended source flux cor¬ 
rection could be significant due to the assumptions that 
the high-J CO transitions are distributed in the same 
way as the the low-J CO lines. If high-J CO transitions 
are only coming from a centralized compact region, 
we are overestimating their flux with our beam correc¬ 
tion method. For this reason, we apply an additional 
15% error to the extended galaxies. There are three 
targets in the sample that have multiple pointings; Arp 
299, NGC 1365, and NGC 2146. In the case of Arp 
299, we use only the pointing for Arp 299 A. For NGC 
1365 we take the average of the northeast and south¬ 
west pointings. This is done since the northeast and 
southwest pointings have approximately a 50% over¬ 
lap in field of view at the center of the galaxy. This 
overlap region is the center of the PACS observations, 
so for comparison, it is best to average the northeast 
and southwest pointings. For NGC 2146, we use only 
the nuclear pointing. 

CO and [Cl] line fluxes were extracted using ver¬ 
sion 1.92 of FTFittei^ a program specifically created 
to extract line fluxes from Fourier transform spectro¬ 
graphs, and are listed in Table [2] This is an interac¬ 
tive data language (IDL) based graphical user inter¬ 
face, that allows the user to fit lines, choose line pro¬ 
files, fix any line parameter, and extract the flux. We 
define a third order polynomial baseline to fit the con¬ 
tinuum for the SLW and SSW separately and derive 
the integrated line intensities from baseline subtracted 
spectra with a simultaneous line fit of all CO, [Cl], 
[Nil] and other bright lines in the spectrum. We use 
a Gaussian line profile, which is based on the assumed 
intrinsic line shape with a width derived from spec¬ 
trally resolved CO 1-0, convolved with the instrumen¬ 
tal line shape, which is a Sine profile. We adopt an er¬ 
ror of 16% for the non-extended galaxy fluxes, which 
encompasses our dominant sources of error, 10% for 
the flux extraction and baseline definition and 6% for 
the absolute calibration uncertainty for staring-mode 
SPIRE FTS observations ( jSwinyard et al.|20l4] l. For 
the case of extended sources, we adopt as already men- 


2 https://www.uleth.ca/phy/naylor/index.php?page=ftfitter 
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tioned, an additional 15% error from the beam size 
corrections, resulting in a total error 30% for the 10 
extended sources. 
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Table 1 

Sample Properties 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

NGC 34 

(IRAS 00085-1223) 

11.49 

0.78 

0.01962 

84.1 

330 

SB 

[OI] 63 . [OI],45.[CII] 
194 - 671 //m 

00* 1 l'"06.67 s -12°06’26.13” 
00*ll'"06.53 s -12°06’24.90” 

1342199416 

1342199253 

KPOT_p vanderw_ 1 
KPOT_p vanderw_ 1 

MCG+12-02-001 
(IRAS 00506+7248) 

11.50 

1.07 

0.01570 

69.8 

200 

SB 

[OI] 63 . [OI] 14 5.[CII] 
194 - 671 //m 

00*54"'03.33 s +73°04'59.83” 
00*54"'03.56 s +73°05'10.38” 

1342193211 

1342213377 

KPOT_p vanderw_ 1 
KPOT_p vanderw_ 1 

IC 1623 

(IRAS 01053-1746) 

11.71 

1.14 

0.02007 

85.5 

250 

sb,agn 

[OI] 63 . [OI] 145 ,[CII] 
194 - 671 //in 

01*07"’46.59 s -17°30’26.46” 
01*07 m 46.74 s -17°30’26.05” 

1342212532 

1342212314 

KPOT_p vanderw_ 1 
KPOT_p vanderw_ 1 

NGC 1068 
(IRAS 02401-0013) 

11.40 

9.07 

0.003793 

15.9 

300 

AGN,SB 

[OI] 63 

[OI] 145 ,[CII] 

194 - 671 //m 

02*42™40.78’ -00°00’47.16” 
02*42™40.73 s -00°00’42.24” 
02*42'"40.92 s -00°00’46.65” 

1342191153 

1342191154 

1342213445 

KPGT_esturm_l K 
KPGT_esturm_l K 
KPGT_c wilsoO 1 _ 1 

NGC 1365 
(IRAS 03317-3618) 

11.00 

4.32 

0.00546 

17.9 

250 

Syl,SB 

[OI] 63 

[OI] 145 ,[CII] 

194 - 671 //m 

03*33"'36.31 s -36°08’16.61” 
03*33 m 36.26 s -36°08’24.33” 
03/133"/36 48 s -36°08’ 19.32” 

1342191295 

1342191294 

1342204020 

KPGT_esturm_lK 
KPGT_esturm_l K 
KPOT_p vanderw_ 1 

NGC 1614 
(IRAS 04315-0840) 

11.65 

1.50 

0.01594 

67.8 

220 

SB 

[OI] 63 . [OI] 145 ,[CII] 
194 - 671 //in 

04*33'"59.79 s -08°34’44.19” 
04*33'"59.85 s -08°34’44.15” 

1342190367 

1342192831 

KPOT_p vanderw_ 1 
KPOT_p vanderw_ 1 

IRAS F05189-2524 

12.16 

0.60 

0.04256 

187 

300 

QSO 

[OI] 63 

[OI] 145 ,[CII] 

194 - 671 //in 

194 - 671 //m 

05*21'"01.24 s -25°21’43.16” 
05*21'"01.28 s -25°21’42.15” 
05*21'"01.42 s -25°21’45.47” 
05*21'"01.42 s -25°21’45.48” 

1342219441 

1342219442 

1342192832 3 

1342192833 3 

KPGT_esturm_l K 
KPGT.esturm.l K 
KPOT.p vanderw. 1 
KPOT_p vanderw_ 1 

NGC 2146 
(IRAS 06106+7822) 

11.12 

6.97 

0.00298 

17.5 

250 

SB 

[OI] 63 . [OI]i 45 ,[CII] 
194 - 671 //m 

06* 18” 35.53 s +78°21’25.39” 
06* 18” 38.07 s +78°21’25.06” 

1342193210 

1342204025 

KPOT.p vanderw. 1 
KPOT_p vanderw_ 1 

NGC 2623 
(IRAS 08354+2555) 

11.60 

1.15 

0.01851 

84.1 

400 

sb,agn 

[Oi] 63 . [OIJ145,[cn] 

194 - 671 //m 

08*38"'24.29 s +25°45'16.72” 
08*38"'24.14 s +25°45'17.34” 

1342208904 

1342219553 

KPOT.p vanderw. 1 
KPOT_p vanderw. 1 

NGC 3256 
(IRAS 10257-4338) 

11.64 

4.61 

0.00935 

38.9 

230 

SB 

[OI] 63 

[OI] 145 ,[CII] 

194 - 671 //m 

10*27"'51.61 s -43°54’15.39” 
10*27"'51.45 s -43°54’21.87” 
10*27"'51.49 s -43°54’16.00” 

1342210383 

1342210384 

1342201201 

KPGT_esturm_lK 
KPGT_esturm_l K 
KPOT_p vanderw. 1 

Arp 299 A 

IC 694 

(IRAS 11257+5850) 

11.93 

4.84 

0.01030 

50.7 

325 

sb,agn 

[OI] 63 

[OI] 145 

[CII] 

194 - 671 //m 

ll*28"'33.41 s +58°33’46.04” 
ll*28"'33.41 s +58°33'46.04” 
ll*28'"33.41 s +58°33'46.04” 
ll*28'"33.41 s +58°33’46.04” 

1342199421 

1342232602 

1342208906 

1342199248 

KPGT_esturm_l K 
OTl_shaileyd_l 
KPGT_esturm_l K 
KPOT.p vanderw. 1 

ESO 320-G030 
(IRAS 11506-3851) 

11.17 

1.67 

0.01078 

41.2 

350 

SB 

[OI] 63 . [OI] 14 5,[CII] 
194 - 671 //m 

11*53'" 11.75 s -39°07’51.75” 
11*53'" 11.52 s -39°07’50.24” 

1342212227 

1342210861 

KPOT.p vanderw. 1 
KPOT.p vanderw. 1 

NGC 4418 
(IRAS 12243-0036) 

11.19 

1.74 

0.007268 

36.5 

163 

Sy2 

[OI] 63 

[oi] 145 ,[cn] 

194-671 //m 

12*26'"54.51 s -00°52’40.77” 
12*26'"54.57 s -00°52’36.93” 
12*26'"54.60 s -00°52’36.54” 

1342187780 

1342210830 

1342210848 

KPGT.esturm.l K 
KPGT.esturm.l K 
KPGT.esturm.l K 

Mrk 231 

(IRAS 12540+5708) 

12.57 

1.48 

0.04217 

192 

200 

QSO 

[OI] 63 

[oi]i4 5 ,[cn] 
194-671 //m 

12*56"' 14.65 s +56°52'24.13” 
12*56"'14.29 s +56°52’23.40” 
12*56"'14.29 s +56°52’26.73” 

1342189280 

1342186811 

1342210493 

KPGT.esturm.l K 
SDP_esturm_3 

KPOT.p vanderw. 1 

IRAS 13120-5453 

12.32 

1.94 

0.03076 

144 

400 

Sy2,SB 

[OI] 63 

[OI] 145,[CII] 

194-671 //m 

13*15'"06.28 s -55°09’24.46” 
13* 15'"06.17 s -55°09’25.38” 
13*15'"06.1 I s -55°09’23.21” 

1342214628 

1342214629 

1342212342 

KPGT.esturm.l K 
KPGT.esturm.l K 
KPOT.p vanderw. 1 

Arp 193 

(IRAS 13183+3423) 

11.73 

8.19 

0.02330 

110 

400 

SB,L 

[OI] 63 . [OI]l4 5 ,[CII] 
194 - 671//m 

13*20"'35.20 s +34°08'24.58” 
13*20"'35.35 s +34°08'23.46” 

1342197801 

1342209853 

KPOT.p vanderw. 1 
KPOT.p vanderw. 1 

NGC 5135 
(IRAS 13229-2934) 

11.30 

0.91 

0.01369 

60.9 

150 

Sy2,SB 

[OI] 63 . [OI]l45,[CII] 
194 - 671 //m 

13*25"'43.96 s -29°50’01.74” 
13‘25'"43.91 s -29°50’00.27” 

1342190371 

1342212344 

KPOT.p vanderw. 1 
KPOT.p vanderw. 1 

ESO 173-GO 15 
(IRAS 13242-5713) 

11.38 

3.61 

0.00974 

34 

200 

SB 

[OI] 63 ,[OI] 145 ,[CII] 
194-671 //m 

13‘27'"24.00 s -57°29’23.63” 
13*27'"23.95 s -57°29’22.89” 

1342190368 

1342202268 

KPOT.p vanderw. 1 
KPOT.p vanderw. 1 

Mrk 273 

(IRAS 13428+5608) 

12.21 

1.05 

0.03736 

173 

520 

SB,Sy2 

[OI] 63 

[OI],45.[CII] 

194 - 671 //m 

13*44'" 42.09 s +55°53'09.14” 
13*44’"41.82 s +55°53'08.75” 
13*44'"42.10 s +55°53'10.50” 

1342207801 

1342207802 

1342209850 

KPGT.esturm.l K 
KPGT.esturm.l K 
KPOT.p vanderw. 1 

Zw 049.057 

CGCG 049-057 
(IRAS 15107+0724) 

11.35 

1.05 

0.01300 

65.4 

200 

SB 

[OI] 63 . [OI]l 4 5,[CII] 
194-671 //m 

15* 13"'13.18 s +07°13'30.71” 
15* 13"'13.10 s +07°13’29.19” 

1342190374 

1342212346 

KPOT.p vanderw. 1 
KPOT.p vanderw. 1 

Arp 220 

(IRAS 15327+2340) 

12.28 

4.87 

0.01813 

77 

504 

sb,agn 

[OI] 63 , [OI]i 45 
[CII] 

194 - 671 //m 

15*34"'57.22 s +23°30’11.06” 
15*34"‘57.21 s +23°30'10.13” 
15*34"'57.11 s +23°30’11.26” 

1342191304 

1342191306 

1342190674 

KPGT.esturm.l K 
KPGT.esturm.l K 
KPGT.c wilsoO 1 _ 1 

NGC 6240 
(IRAS 16504+0228) 

11.93 

1.10 

0.02448 

116 

500 

sb,agn 

[OI] 63 

[OI] 145,[CII] 

194-671 //m 

16*52"'59.10 s +02°24'03.58” 
16*52'"59.10 s +02°24'02.79” 
16*52'"59.01 s +02°24’03.27” 

1342216622 

1342216623 

1342214831 

KPGT.esturm.l K 
KPGT.esturm.l K 
KPOT.p vanderw. 1 

IRAS FI7207-0014 

12.46 

1.56 

0.04281 

198 

620 

SB,L 

[OI] 63 

17*23'»2i.83 s -00°16’59.82” 

1342229692 

KPGT.esturm.l K 
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2.3. Herschel/PACS observations 


We have obtained observations of the [OI] 63 pm 
([OI] 63 ), [OI] 145 pm ([OI] 145 ), and [CII] 157 /urn 
emission lines with the Integral Field Spectrometer 
of the Photodetector Array Camera and Spectrometer 
(PACS, Poglitsch et al.|2010| > on board Herschel Space 
Observatory for every object in the HerCULES sam¬ 
ple. The data presented here were obtained as part of 
the Herschel program KPOT pvandcrw 1 (PI: P van 
der Werf), complemented by observations from other 
programs. The observations and program IDs of the 
[CII] and [OI] lines are listed in Table [T] 


The data were downloaded from the Herschel Sci¬ 
ence Archive and processed using HIPE vll.O. Stan¬ 
dard processing steps including timeline deglitching, 
application of the Relative Spectral Response Func¬ 
tion and detector flat fielding, and subtraction of the 
on and off chop positions, gridding along the spec¬ 
tral axis, and combination of the nod positions. With 
the exception of Arp 299, the objects are all centered 
on the 974 central spaxel of the 5 by 5 PACS array, 
observed in staring mode. The fluxes are extracted 
from the central spaxel, using the extractSpaxelSpec- 
trum routine, and referenced to a point source. We 
use the pointSourceLossCorrection routine to capture 
any additional flux that may not be captured in the 
central spaxel. Finally, version 3.10 of SPLAT as 
part of the STARLINK software package (http://star- 
www.dur.ac.uk/~pdraper/splat/splat-vo/) was used to 
subtract the baseline from each observation, and iso¬ 
late the desired lines, in the case of PACS range spec¬ 
troscopy. The reduction steps were the same for both 
the PACS range and line spectroscopy, two different 
observing modes of PACS. 


Arp 299 was observed in the mapping mode and re¬ 
duced using the standard pipeline reduction. The inte¬ 
grated flux for Arp 299 A (presented in this paper) was 
calculated by summing the flux within a 25” aperture 
centered on Arp 299 A SPIRE pointing. 

In order to extract the line parameters from the 
PACS observations, we first integrate over the baseline 
subtracted spectrum and then we fit a Gaussian pro¬ 
file to the baseline subtracted flux. In some sources, 
the [OI ]63 line shows a double-peaked profile, where 
the flux at the central wavelength is diminished, which 
could indicate Keplerian rotation. However, if this 
were the case, then we would expect a similar pro¬ 
file in the [OI] 145 pm line and possibly the other fine 
structure lines as well, which is not seen. The spec¬ 


tral resolution of PACS at 145 pm is more than suffi¬ 
cient to resolve the ~ 0.2 pm separation between the 
two peaks in the [OT ] f ,3 profile (Figure |T]). Therefore, 
we conclude that this double-peaked profile is due to 
absorption in the center of the profile by colder fore¬ 
ground gas. We note that [OI]® absorption is due to 
O in the ground state while absorption at 145/rm re¬ 
quires O to be at a state having an energy of 226 K 
above the ground state. Therefore, in cool or moder¬ 
ate density gas the [OI]i 45 line will not show an ab¬ 
sorption feature, even if the [OI ]63 line does. This 
same effect has been noted in Arp 220, which shows 


the [OI]® in full absorption (Gonzalez-Alfonso et al. 


[20l2 ). In the case of NGC 4418 and Zw 049.057, the 
[OI ]63 line has an inverse P Cygni profile, suggesting 
that the absorbing foreground gas is flowing into the 
nuclear region. The three example galaxies for which 
the spectra are shown in Figures 1-3 display increas¬ 
ing absorption of the [OI ]63 pm line. In NGC 7552, 
a face-on starburst galaxy, the profile remains Gaus¬ 
sian, while in Mrk 331, a late-stage merger, there is a 
strong dip in the middle of the profile. IRAS F17207- 
0014 is known for being one of the coolest ULIRGs, 
here absorption dominates the [OI ]63 emission. For the 
[OI]® profiles that show an absorption feature, we fit 
the Gaussian only to the wings of the emission pro¬ 
file and state the flux in parentheses. The Gaussian- 
fit flux is only valid if the true line profile is Gaus¬ 
sian. We suggest this is a more robust estimate of the 
true integrated flux of the [OI]® line emerging from 
the warm nuclear region, since in many cases, the ab¬ 
sorption dominates the profile. We note that using a 
Gaussian profile to extrapolate the line flux requires 
an assumption that the location of the line center (a 
free parameter) is in the middle of the profile, which 
may not be accurate, especially in the case of IRAS 
F17207-0014, or any other galaxies with an asymmet¬ 
ric line profile. We have tested the relations presented 
in the rest of this paper with both the integrated flux 
and the Gaussian fit, and find it does not strongly af¬ 
fect the results. Both the observed line fluxes and the 
gaussian-fit line fluxes, stated in the parentheses, are 
also presented in Table [2] 


3. Results 

3.1. Spectra and line fluxes 

All SPIRE CO and [Cl] line fluxes are listed in Ta¬ 
ble^ We present three examples of galaxy spectra ob¬ 
tained with SPIRE in the top panels of Figure |T][2j and 
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Table 1— Continued 


1 

2 

3 

4 

5 

6 

7 

8 


9 

10 

11 








[OI] 145 .[CII] 

17*23' 

"31.84 s -00° 16’57.60” 

1342229693 

KPGT_esturm_l K 








194 - 671 /im 

17*23' 

"21.93 s -00°17’01.10” 

1342192829 

KPOT_p vanderw_ 1 

IRAS FI8293-3413 

11.88 

1.82 

0.01817 

86 

270 


[OI] 63 , [OI] 145 ,[CII] 

18*32' 

"41.34 s -34°IT36.90” 

1342192112 

KPOT_p vanderw _ 








194 - 671 /im 

18*32' 

"41.17 s -34° 11’27.23” 

1342192830 

KPOT_p vanderw_ 1 

IC 4687/6 

11.62 

0.84 

0.01735 

81.9 

230 

SB 

[OI] 6 3 

18*13' 

"39.94 s -57°43’49.66” 

1342239740 

OTIJarmus.l 

(IRAS 18093-5744) 







[CII] 

18*13' 

"39.80 s -57°43’35.71” 

1342239739 

OT1 Jarmus.l 








194 - 671 /im 

18*13' 

"39.50 s -57°43’31.05” 

1342192993 

KPOT_p vanderw. 1 

NGC 7469 

11.65 

1.32 

0.01632 

70.8 

300 

Syl,SB 

[OI] 6 3 

23*03" 

'15.47 s +08°52’37.05” 

1342187847 

KPGT_esturm_l K 

(IRAS 23007+0836) 







[OI]i 45 ,[CII] 

23*03" 

'15.83 s +08°52’28.52” 

1342211171 

KPGT_esturm_l K 








194 - 671 //m 

23*03” 

'15.87 s +08°52’28.15” 

1342199252 

KPOT_pvanderw_ 1 

NGC 7552 

11.11 

3.64 

0.00537 

23.5 

180 

SB 

[OI] 63 

23*16' 

"10.10 s -42°34’53.89” 

1342210400 

KPGT_esturm_l K 

(IRAS 23134-4251) 







[OI] 145 ,[CII] 

23*16' 

"10.59 s -42°35’05.73” 

1342210399 

KPGT_esturm_l K 








194 - 671 /im 

23*16' 

"10.73 s -42°35’06.02” 

1342198428 

KPOT_p vanderw_ 1 

NGC 7771 

11.40 

1.08 

0.01427 

61.2 

250 

SB 

[OI] 63 , [OI] 145 ,[CII] 

23*51" 

'24.83 s +20°06’42.33” 

1342197839 

KPOT_pvanderw_ 1 

(IRAS 23488+1949) 







194 - 671 /im 

23*51” 

'24.72 s +20°06’42.11” 

1342212317 

KPOT_p vanderw_ 1 

Mrk 331 

11.50 

0.87 

0.01790 

79.3 

215 

SB 

[OI] 63 , [OI] 145 ,[CII] 

23*51" 

'26.76 s +20°35’09.83” 

1342197840 

KPOT_p vanderw. 1 

(IRAS 23488+2018) 







194 - 671 /im 

23*51" 

'26.65 s +20°35’ 10.42” 

1342212316 

KPOT_pvanderw_ 1 


Note.— 

Column 1: Object name. 

Column 2: log(LiR/L 0 ) from |Armus et al.|{2009f . Observations use the cosmological param eters Ho=70 km s _I Mpc, £l vaccum =§.12, and Gl matter =0.2&. 

Column 3: Far infrared flux (FIR) calculated using the IRAS definition jHelou et al.|l985| in 10 -12 W m -2 . 

Column 4: Redshift z from NED. 

Column 5: Luminosity distance in Mpc from |Armus et al.|j2009) . 

Column 6: CO 1-0 full width to half power line wI3th in km s . 

Column 7: Galaxy classification from NED SB=Starburst, L=LINER, AGN=Active Galaxy Nucleus, Syl=Seyfert 1, Sy2=Seyfert 2, QSO=Quasi-Stellar Object. 
Column 8: Line names. 

Column 9: Pointing Coordinates. 

Column 10: Observation ID (OBSID). 

Column 11: Program ID. 

a The two SPIRE observations of IRAS F05189-2524 were combined using HIPE average (avg) task. 
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[3] for NGC 7552, Mrk 331, and IRAS F17207-0014, 
respectively. It is important to note that the baseline 
ripple seen in the SPIRE FTS spectra is due to both 
the sine profile of the strong CO transitions and the 
noise. Since in this paper we only discuss the neutral 
gas cooling, we do not present fluxes of [Nil] (which 
originates in ionized gas) or the molecular lines other 
than CO, which are irrelevant to the total neutral gas 
cooling. A comprehensive set of fluxes will be pre¬ 
sented in Van der Werf et al. (in prep). In addition, 
the HerMES team is planning to publish a formal data 
paper using HIPE vl2.0 with a detailed error analysis 
in the near future. 

In the bottom row of the spectra in Figures |T| [2] 
and [3j the PACS line profiles of the three sample 
galaxies are presented (NGC 7552, Mrk 331, IRAS 
F17207-0014). 
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Rest Wavelength[|im] 



400 600 800 1000 1200 1400 1600 

Observed Frequency [GHz] 


[OI] 63 [01], 45 [Cl l]l 57 



Wavelength [/xm] Wavelength [/j,m] Wavelength [^m] 


Fig. 1.— Top: Herschel SPIRE spectrum for NGC 7552 in observed frequency. CO lines are indicated in blue and 
tine structure lines in pink. The vertical black dashed line near 920 GHz separates the SLW and SSW arrays. The 
inset shows a magnified version of the CO J=7-6 and [CI] 37 o transitions. Bottom: Baseline subtracted Herschel PACS 
observations of [OI]g 3 , [OI] 145 , and [CII] presented in black, with superimposed Gaussian line fits shown in red. 
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Rest Wavelength[pm] 



400 600 800 1000 1200 1400 1600 

Observed Frequency [GHz] 

[OI] 63 [OI] 145 [Cll] 157 



Wavelength [/^m] Wavelength [/xm] Wavelength [/im] 


Fig. 2.— Top: Herschel SPIRE spectrum for Mrk 331 in observed frequency. CO lines are indicated in blue and 
fine structure lines in pink. The vertical black dashed line near 920 GHz separates the SLW and SSW arrays. The 
inset shows a magnified version of the CO .1=7-6 and [CI] 37 o transitions. Bottom: Baseline subtracted Herschel PACS 
observations of [OI] 63 , [OI] 145 , and [CII] presented in black, with superimposed Gaussian line fits shown in red. 
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Rest Wavelength[)im] 

700 600 500 400 300 ' 200 



Observed Frequency [GFIz] 

[OI]63 [OI]l45 [Cl l]l 57 



Wavelength [/j,m] Wavelength [nm\ Wavelength [^,m] 


Fig. 3.— Top: Herschel SPIRE spectrum for IRAS F17207-0014 in observed frequency. CO lines are indicated in 
blue and fine structure lines in pink. The vertical black dashed line near 920 GHz separates the SLW and SSW arrays. 
The inset shows a magnified version of the CO .1=7-6 and [CI] 37 o transitions. Bottom: Baseline subtracted Herschel 
PACS observations of [OI] 63 , [OI] 145 , and [CII] presented in black, with superimposed Gaussian line fits shown in red. 
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Object 

CO J= 1 -0 

Beam (”) 

Ref. 

CO J=2-l 

Beam (”) 

Ref. 

CO J=3-2 

Beam (”) 

Ref. 

NGC 34 

0.74 

48” 

H98 

3.06 

23” 

H98 

- 



MCG+12-02-001 

- 



- 



- 



IC 1623 

2.60 

22 ” 

P12 

17.11 

52” 

104 

37.5 

ii” 

P12 

NGC 1068 

10.82 

22 ” 

P12 

86.37 

14” 

P12 

196.05 

ii” 

P12 

NGC 1365 

13.16 

55” 

G04 

18.70 

23” 

H98 

53.90 

14” 

114 

NGC 1614 

0.79 

21 ” 

G04 

1.54 

22 ” 

K13 

- 



IRAS F05189-2524 

0.18 

22 ” 

P12 

0.96 

14” 

P12 

2.83 

11 ” 

P12 

NGC 2146 

9.53 

21 ” 

G04 

6.91 

12.5” 

B93 

84.46 

21 ” 

M99 

NGC 2623 

0.61 

22 ” 

P12 

2.01 

14” 

P12 

6.94 

11 ” 

P12 

NGC 3256 

3.30 

44” 

A95 

53.00 

22 ” 

A95 

- 



Arp 299 

2.23 

22 ” 

P12 

- 



49.66 

11 ” 

P12 

ESO 320-G030 

0.68 

48” 

M90 

- 



- 



NGC 4418 

0.50 

22 ” 

P12 

- 



11.37 

11 ” 

P12 

Mrk 231 

0.32 

22 ” 

P12 

2.32 

14” 

P12 

6.27 

11 ” 

P12 

IRAS13120-5453 

- 



- 



- 



Arp 193 

0.73 

22 ” 

P12 

6.37 

14” 

P12 

13.55 

11 ” 

P12 

NGC 5135 

1.45 

22 ” 

P12 

9.35 

14” 

P12 

22.25 

11 ” 

P12 

ESO 173-GO 15 

- 



- 



- 



Mrk 273 

0.30 

22 ” 

P12 

2.00 

14” 

P12 

5.35 

11 ” 

P12 

Zw 049.057 

0.45 

22 ” 

P12 

4.58 

14” 

P12 

8.06 

11 ” 

P12 

Arp 220 

1.58 

22 ” 

P12 

8.49 

14” 

P12 

41.53 

11 ” 

P12 

NGC 6240 

1.21 

22 ” 

P12 

11.17 

14” 

P12 

36.00 

11 ” 

P12 

IRAS F17207-0014 

0.59 

22 ” 

P12 

5.06 

14” 

P12 

13.22 

11 ” 

P12 

IRAS FI 8293-3413 

2.23 

55” 

G04 

- 



- 



IC 4687/6 

0.42 

48” 

A07 

2.65 

23” 

A07 

- 



NGC 7469 

1.12 

22 ” 

P12 

6.72 

14” 

P12 

18.12 

11 ” 

P12 

NGC 7552 

3.10 

48” 

C92 

21.00 

22 ” 

A95 

26.00 

15” 

114 

NGC 7771 

1.33 

55” 

S91 

- 



10.12 

23” 

N05 

Mrk 331 

1.27 

55” 

G04 

- 



- 




Table 3: Ground based CO integrated fluxes from the litera¬ 
ture. All units are in 10 18 W m~ 2 . References are as follows: 
A95 =|Aalto et al.|(l995bl , A07=|Albrccht et al .1(2007) . G04 


[&So 

(2004 


omon M2004| , H89= |Heckman et al. 1(1989) . I04=[i ono et al. 


et al.|i 


1999 


2012 


I14=Israel, F.R, 2014 A&A to be submitted, K13=|Kbnig 


20131 , M90=|Mirabel et ah j 19901 , M99=jMauersberger et al. 
i, N05= |Narayanan et al.| ~ 20051, P12= [Papadopoulos et al. 


, and references therein, S91=|Konig et al ,|{20 1 3}. 
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3.2. Classification of CO ladders 

In Figure [4] we present the CO ladders of the full 
HerCULES sample. We have collected the available 
ground based observations of CO J=l-0, 2-1, and 3- 
2, whose fluxes and references are listed in Table [3] 
Where necessary we have converted the ground-based 
measurements to the cosmology adopted here (Section 
1). In order to compare these CO ladders directly, we 
have normalized the ladders by the integrated CO flux 
summed from .1=4-3 through 1=13-12, to focus on the 
relative behavior of the higher-J transitions since we 
do not have CO J=l-0 data for all sources. For galax¬ 
ies with non detections for any CO transitions, we use 
the linearly interpolated value. The CO ladders are 
separated into three classes based on the parameter a , 
where: 


CO emission, when also considering additional con- 


straints. Arp 220 fRangwala et al. 

2011), Arp 299 

(Rosenberg et al. 2014b 

), NGC 253 ( 

Rosenberg et al. 

2014a), and NGC 6240 

Meijerink et al. 20131 require 


mechanical heating to reproduce the high-J CO lines, 
while Mrk 231 ( van der Werf et al.|2010b| > and NGC 


1068 (jSpinoglio et al.|2012 1 require X-rays to directly 


heat the gas in order to reproduce the observed molec¬ 
ular emission. This trend suggests that when dealing 
with highly excited CO ladders, such as in Class II 
and especially Class III objects, there is an additional 
heating mechanism necessary to explain the observed 
molecular emission. We will explore this issue for the 
full sample in the next section. 


4. Analysis 


LcOj= 11-10 +TcOj = 12-11 + L C Oj= 13-12 

a - ---. ( 1 ) 

LcOj= 5-4 + L C Oj= 6-5 + LcOj= 7-6 

and Leo — ^xD^Fco, with F co in [W m ‘ 2 ] (TableJ^J) 
and luminosity distance, D~ in [m] (listed in Table fill. 
Here we use three transitions of both the mid- and 
high-J CO transitions to help prevent noise or a non¬ 
detection of one of these lines from dominating a. We 
define the three classes as: 


• Class I: a < 0.33 

• Class II: 0.33 < a < 0.66 

• Class HI: a > 0.66 


The definition of the classes is quantitatively arbitrary, 
but chosen to reflect similarities in the spectral line en¬ 
ergy distributions, which is illustrated in Figure [4] In 
the case that we do not have any observations of the 
low-J transitions, we do not plot any low-J fluxes. The 
parameter a is based on the ratio of three high-J CO 
lines to three mid-J CO lines, which essentially defines 
the drop-off slope of the CO ladder from J=5-4. Thus, 
the steepest drop-offs are in Class I, while the flattest 
ladders are in Class III. Class II consists of objects that 
peak around J=6-5, but do not fall off as steeply as 
those of Class I. Our three example galaxies were se¬ 
lected to fit into these categories, with NGC 7552 as 
a Class I, Mrk 331 as a Class II, and IRAS F170207- 
0014 as a Class III object. 

We note that the CO ladders for many Class II and 
III (the excited classes) objects have been published. 
In the case of all of these sources, heating mechanisms 
besides UV heating are required to explain the high-J 


We will combine our PACS and SPIRE observa¬ 
tions of all the major neutral gas cooling lines ([OI], 
[Cl], [CII], and CO), with ancillary data obtained as 
part of The Great Observatories All-Sky LIRG Sur¬ 
vey (GOALS) (Armus et a H2009| l, including the ma¬ 
jor PDR coolant [SiH] at 34.8pm. We choose to in¬ 
clude [Sill] since this is a known PDR coolant that 
has strength on order of [CII], thus is an important el¬ 
ement of the cooling budget. We will not be dealing 
with any ionized gas coolants such as [Nil], [OIII]. We 
use the IRAS definition of the far infrared flux (FIR) 
as FIR= 1.26 x 1(T 14 (2.58 S 60 ^+Sioo Mm )[W m~ 2 ] , 
where S v is in units of Jansky [Jy] ([Helou et al. |l985ji. 
We then use the luminosity distances (D/J, from |Ar-| 
|mus et ah] ( |2009j ). to define the far infrared luminos¬ 
ity (TfirX making our calculations directly compara¬ 
ble to |Stierwalt et al.| ( |2013[ ). When we refer to the CO 
flux, we use the sum of the line fluxes from CO J=4-3 
through J= 13-12. 


4.1. Warm gas tracers 

In the local universe, the ratio of IRAS 60/100 pm 
flux densities correlates with infrared luminosity (e.g., 
Chap man et al.||2003| i. It is of interest to determine 
whether the IRAS 60/100 pm ratio (a proxy for dust 
temperature) or Lfir (dust luminosity) correlates better 
with the degree of CO excitation, as parametrized by a. 
Since a represents a proxy for the slope of the CO lad¬ 
der above J=5, it traces the relative brightness of high- 
J lines in comparison to mid-J lines, allowing a rough 
estimate of overall CO excitation. When a is small, 
the CO excitation is low and the CO SLED is highly 
peaked and when a is large, the CO SLED is flattened 
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Fig. 4.— CO spectral line energy distributions for the full HerCULES sample divided into three classes. Class I (left 
panel) includes galaxies with a < 0.33, Class II (center panel) is where 0.33 < a < 0.66, and Class III (right panel) 
is where 0.66 > a; a is defined in Eq. |T| in Section 3.2 Where we do not have data, we have linearly interpolated 
between the neighboring transitions. If we lack all three ground based transitions, we do not plot any low-J fluxes. 
The ground based transitions have been plotted in gray, to emphasize the Herschel line fluxes. 


and the excitation is high, indicating significant emis¬ 
sion by warm and dense molecular gas. We compare 
the Sgo^m/Sioo/jm ratio from Sanders et al. (2003) and 
the Lfir to the molecular gas excitation (a). In the 
left panel of Figure [5] the excitation (a) is plotted as a 
function of S^m/S lOO/m, where each square point is a 
galaxy in our sample. The best-fit power law is shown 
with a red dashed line. In the right panel, a similar cor¬ 
relation between Lfir and a is shown, with the best-fit 
power law plotted with a red dashed line. We see that 
although both panels show a positive trend, the corre¬ 
lation found with the S^O/jm/Swo/jm ratio is tighter than 
that seen with the Lfir. The molecular gas excitation 
to infrared color relation has a correlation coefficient 
of an r=0.81, while the excitation to Lfir relation has 
a correlation coefficient of an r=0.52. Although the 
excitation to infrared color relationship is significantly 
correlated, we find three outliers, based on the largest 
Euclidean distance and shown in red in the left panel 
of Figure [5j, Arp 299, NGC 1614, and NGC 2623. We 
also select the three farthest outliers, also based on the 
Euclidean distance, in the a to L F i R relation (Figure [5] 


Arp 299, NGC 4418 , and IRAS FI8293-3413), also 
plotted as red points in the right panel. We use the tra¬ 
ditional definition of Euclidean distance: \\u - v|| 2 . To 
test the strength of the correlations, we refit a power 
law excluding the three outliers in each plot, which are 
marked in red. We remove the outliers to test if the 
correlation between the CO excitation and the IRAS 
colors is still higher than that of the L/./«, or if it is just 
the outliers affecting the coefficient. The new best-fit is 


plotted as the blue solid line in Figure [5] Although the 
exclusion of these points does not result in a significant 
change in the best fit in either case, it does improve 
the correlation coefficients, resulting in a correlation 
coefficient of an r=0.88 for the molecular gas excita¬ 
tion to infrared color correlation, and an r-0.10 for the 
molecular gas excitation to Lfir relation. Physically, 
this suggests that the presence of warm, dense molecu¬ 
lar gas, is correlated with the presence of warm dust. It 
is, however, important to note that, once removing the 
outliers, the correlation between a and L F ir is also sig¬ 
nificant, with a 1.6% probability of this relation being 
spurious before removing the outliers. Further, since 
the IRAS 60/100 pm ratio is shown to correlate with 
Lfir, these two quantities are likely related by under¬ 
lying variables, making this correlation difficult to in¬ 
terpret. We note however that Lu et al. ( 2014) 1 com¬ 
pared specific CO line transitions normalized by FIR 
and the IRAS 60/100 /im flux ratio. They find that as 
the CO gas becomes warmer, the 60/100 /im ratio also 
increases, which is in good agreement with our results. 


4.2. Cooling budget 

We can calculate the neutral gas cooling budget 
in each galaxy by summing the luminosities of the 
[OI], [Cl], [Sill], [CD], and CO lines, since these are 
the main neutral gas coolants in the mid- and far- 


infrared regime. We take the [Sill] fluxes from Inami 


|et al.| ( |2013| ), which were observed with the Spitzer 
IRS instrument in the long wavelength, high resolu- 


16 



























Fig. 5.— Left panel: The gas excitation (a, for a definition see Section [L2| plotted against the IRAS infrared color 
(Sf,opm/S [oo/rmf each square is a galaxy in our sample. The red dashed line and text represents the least squares fit 
and correlation coefficient (r) for the full sample. The three red circles represent the the most extreme outliers (largest 
Euclidean distance value) in the relation (Arp 299, NGC 1614, and NGC 2623). The blue solid line and text represent 
the best fit and correlation coefficient for the sample excluding the three most extreme outliers, also using the largest 
Euclidean distance. Right panel: The same as the left panel but the gas excitation is plotted against the Lfir, calculated 
in Section[4] The three outliers in this case are Arp 299, NGC 4418, and IRAS F18293-3413. 


tion mode. We note that the slit size for Spitzer IRS- 
LH is ll.l”x22.3”, thus galaxies that are more ex¬ 
tended than 11.1” will be missing some [Sill] flux. 
There are five affected galaxies, Arp 193, Arp 299, 
IRAS 13242-5713, NGC 1365, and NGC 2146. In 
Figure [ 6 j we present the percentage of cooling con¬ 
tributed by each emission line as a function of Lfir. 
The percentage cooling for each species is calculated 
by comparing the luminosity of a species to the to¬ 
tal summed luminosity of the [OI], [Cl], [Sill], [CD], 
and CO lines. We exclude three galaxies based on 
their [OI]g 3 profiles, Arp 220 that is fully in absorp¬ 
tion, and NGC 4418 and Zw 049.057, which are both 
heavily absorbed and show inverse P-Cygni profiles 
(Gonzalez-Alfonso et al.] |2012j >. For each galaxy, the 
percent of cooling contributed by [CI ]609 and [Cl] 370 is 
plotted as a yellow circle, CO J=4-3 through J=13-12 
in red, [Sill] in blue, and the combined cooling of [OI] 
63/7111, [OI] 145pm, and [CII] is plotted in green. In the 
bottom panel of Figure [ 6 ] we separate the cooling con¬ 
tributions of [OI] and [CII]. For the five galaxies that 
are larger than 11.1”, we expect a higher [Sill] contri¬ 
bution than shown in Figure [ 6 ] The solid lines show 
the mean percent-cooling for each emitting species. 


For example, the neutral atomic carbon is responsi¬ 
ble for no more than 2 % of the total cooling with an 
average cooling contribution of 1.5%, while CO con¬ 
tributes a mean of 10.8%, and [Sill] contributes 24.2%. 
There are two galaxies with exceedingly high CO cool¬ 
ing percentages, namely IRAS F05189-2524 and Mrk 
231, the two strongest AGN in the sample. Very high 
CO cooling percentages have also been noted in the 
massive Galactic star forming region W3 with 32% to¬ 
tal gas cooling ( |Kramer et al.||20 04[>, and even higher 
percentages in DR21 (Jakob et al.||2007 White et al. 


2010 1 . In both cases, the high CO percentage is at¬ 


tributed to self absorption of the [OI ]63 line, yet we 
calculate our fluxes both with the Gaussian and ob¬ 
served [OI ]63 line fluxes and see little change. There¬ 
fore, for IRAS F05189-2524 and Mrk 231, we believe 
the high CO cooling percentage is a true effect, and 
not one dependent on the [OI] (,3 absorption. The most 
efficient coolants are [OI] and [CII], which together 
provide a mean of 63.7% of the total gas cooling bud¬ 
get. Separating their cooling contributions, [CII] pro¬ 
vides a mean cooling percentage of 33.6% and [OI] 
cools a mean of 30.1% of the gas. The mean cooling 
percentages and their standard deviations are shown in 
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Table [4] Inspection of Figure [6] shows that the outliers 
are randomly distributed, and there is no clear trend of 
outliers as a function of L F ir. 

Table 4: Mean and standard deviation of percent cool¬ 
ing contribution. 


Line 

Mean 

Std. Dev. 

[CII] + [OI] 63+ 145 

63.7 

14.3 

[CII] 

33.6 

9.1 

[OI]63+i45 

30.1 

11.8 

[Sill] 

24.2 

9.9 

ELcOj 

27=i [ci]j 

10.8 

10.0 

1.5 

0.9 


It is interesting to note how constant each cooling 
range is as a function of Lfir. It is a well known phe¬ 
nomenon that as far infrared luminosity increases, an 
apparent [CII]/FIR ratio decreases, which is what we 
call the [CII] deficit. This is observed in various envi¬ 
ronments in the local Universe ( Malhotra et a l. |2001[ 
|Luhman et al.|[2003 DIaz-Santos et al.||2013 1. In ad¬ 
dition, Gracia-Carpio et al. (2011) find a fine-structure 
line deficit in [Nil] and [OI] as well, such that both the 
[NII]/FIR and [OI]/FIR ratios decrease as a function 
of far infrared luminosity. One corollary to this fine- 
structure line deficit is that other (molecular) coolants 
could become more efficient at higher Lfir. However, 
as shown in Figure [6] the relative efficiencies of each 
coolant remain mostly constant. In order to better un¬ 
derstand this phenomenon, we plot the [CII] deficit 
([CII] to FIR ratio), as a function of Lfir in the top 
left corner of Figure [7] For comparison, we plot the 
same ratios but for the [OI] 63 +i 45 , [Sill], CO, and [Cl] 
lines, where the CO lines encompass all transitions 
4< J upp <13 and [Cl] is the sum of both [CI ]609 and 
[Cl] 370 - We see a trend that is consistent with a deficit 


in [CII], agreeing with the observations from DIaz- 


Santos et al. (20131. Since our sample spans a smaller 


luminosity range than DIaz-Santos et al. ( 2013] >, our 
deficit only spans a factor of 2 in the [CII]/FIR ratio, 
while in the GOALS sample, it spans an order of mag¬ 
nitude. Thus, we plot the deficit in linear space and see 
clearly that the only points with a high [CII]/FIR ratio 
are those of lower Lfir- We see a similar trend in [Sill] 
and [Cl], where the latter has also been observed from 
ground based observatories (Gerin & Phillips 1998 
|2000| ). The [OI] shows a tentative line deficit, which 
becomes more obvious when we include galaxies of 
higher IR luminosities ( Gracia-Carpio et al.|2011 1 . Al¬ 
though we find evidence consistent with line deficits in 


[CII], [OI], [Sill], and [Cl], there is no evidence for a 
deficit in CO, which shows a flat distribution over Lfir, 
further strengthening the results from Lu et al. ( 2014| >. 
In all panels, NGC 6240 is an extreme outlier, as was 
also noted in |Lu et ak] ( |2014| >. 


The difference between the molecular and fine 
structure emission can be understood in terms of heat¬ 
ing mechanisms. The fine structure lines are origi¬ 
nating from regions that are heavily affected by UV 
photons, at the edges of PDRs. As shown in the PDR 
models from |Kaufman et al.| ( |1999] >, as the radiation 
field and density increase, the fine structure line emis¬ 
sion is expected to weaken compared to the far infrared 
flux. However, this does not apply for the molecular 
gas (CO), where we see no line deficit. This result 
may not be surprising since CO traces the molecular 
gas deeper in molecular clouds, where the UV field 
is significantly attenuated. We can test the gas versus 
dust cooling efficiency by plotting the ratio of the to¬ 
tal gas cooling to the far infrared flux as a function of 
Lfir, which is shown in the bottom right panel of Fig¬ 
ure [7] Here, it is clear that the trend is decreasing in a 
very similar manner to that of the fine structure lines, 
which is expected since [CII],[OI], and [Sill] domi¬ 
nate the cooling. Therefore, we can say that the fine 
structure line deficit is actually a gas cooling deficit 
in comparison to the FIR flux that is a result of UV 
heating becoming a more efficient coolant of dust in 
exceedingly extreme environments. It is critical to note 
that the fine structure line fluxes are not decreasing in 
absolute flux, but their relative contribution in com¬ 
parison to the warm dust (measured with the FIR flux) 
is decreasing, due to the warm dust becoming increas¬ 
ingly efficient at a faster rate than the gas coolants. 
This can also be understood as dust being heated more 
efficiently, such that the fine structure lines do not ac¬ 
tually show a deficit, which agrees with the results of 
), where they find that grains 
are ’stealing’ photons from the gas. 


DIaz-Santos et al. ( j2013 


However, if the neutral gas cooling efficiency is de¬ 
creasing, and we observe a tentative line deficit in all 
species except for CO, then we would expect the per¬ 
centage of CO cooling to increase slightly as a function 
of Lfir, to compensate for the decreasing fine structure 
line cooling efficiency. We do see a slight increase in 
the percentage of CO cooling, but the trend is tentative 
and within the general scatter of the other CO cool¬ 
ing percentages. In addition, the galaxies with high 
Lfir have flatter CO ladders, meaning there is non- 
negligible flux in high-I (J> 13) transitions, that we 


18 
















































80 


60" 


40 


CM CD CO 
ID CD 
ID i- O 
N CVJ i- 

o oo 
0 00 
z z z 


-9 

£*• CO 

oo 

00 ) 


V ow cd co 

CM CDi- Lf) CM 

-M-CD CM CDco 
+ COCO CO I''-’- CO CMCM 

oo^Sooo o® 
oggtooo ta, 
s z^y zz z zy 





1 




si 

t 

| 



>< 








i * 

r ; 

V 

r 


r i 

S' 

! 

'i 

★ 

» 

jr 





♦ 

ijr 

i 


r 


• i 


_ 

_ 


• 




i; 




■ 

• 





* ki 

u 

* 

I 


af 

r 

£ 



• : 

• 

n 


11.C 

1 

1.5 

12.0 


=5 20" 
o 
O 

| ot 

Lj 

270 

D) 

C 

1 60 
O 
0 ) 

^ 50 
75 

,2 40 

D) 

= 30 
o 
o 
O 

a> 20 


log(Lnfl) [L 0 ] 




4» ♦♦ ■ 




* 


11.0 


11.5 

iog(i™) [L] 


12.0 


[CIIMOI] 

[Sill] 

Em co, 
EmICI], 


[Cll] 

[OI]e3+145 


Fig. 6 . — Top: Percentage of line cooling of the total gas cooling from each observed species as a function of Lfir- 
Given are absolute, not cumulative percentages. The yellow triangles are the percentage of cooling from [CI ]609 
and [CI] 37 o, the red circles are the cooling from CO (4< J upp < 13), the blue stars are the cooling from [Sill], and 
the green pentagons are from [OI] 63 + 145 +[CII]. The colored lines represent the mean cooling percentages for each 
coolant. We exclude galaxies that show [OI] in absorption (Arp 220), or show a complex line profile (NGC 4418 and 
Zw 049.057). Bottom: Same as top but separating the cooling contributions of [CII] (cyan squares) and [OI] 63 +i 45 
(magenta diamonds). 
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Fig. 7.— The ratio of line fluxes to far infrared flux (FIR) plotted as a function of far infrared luminosity (Lfir) (both 
defined in Section [4]> for [CII] (top left), [OI] (top right), [Sill] (middle left), CO (middle right), [Cl] (bottom left), 
and the total cooling (bottom right). The [OI] flux is derived from the sum of [OI ]63 and [OI]i 45 , the CO flux is the 
sum of line transitions from 1=4-3 to J=13-12, and the [Cl] flux is the sum of [CI ]609 and [Cl] 370 . The total cooling 
refers to the total neutral gas cooling and includes all aforementioned fluxes. We exclude Arp 220, NGC 4418, and 
Zw 049.057 from the [OI] deficit and total neutral gas cooling deficit plots, since we only observe these lines in full or 
partial absorption and can thus not get an accurate flux estimate. 
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do not account for in our cooling budget. This missed 
high-J flux would increase the CO cooling percentage, 
and could make the CO cooling budget increase for 
higher luminosity sources. Since we only have four 
galaxies in our sample within the range that we would 
expect to see an elevated CO cooling percentage, we 
would need to increase our sample size in the high lu¬ 
minosity regime to determine if this trend is real. 

These results have several implications. Firstly, it is 
unlikely that the [CII] line deficit could be caused by 
high dust opacities at 158 /im, since the line deficit is 
also observed in lines of much longer wavelengths, no¬ 
tably the [Cl] line. Similarly, the line deficits are not 
likely caused by the line being optically thick, since 
the deficit is also observed in the optically thin [Cl] 
lines, which agrees with the results of 
et al.j ( j2013| ). A likely explanation is 
tion dominated by ultracompact HII regions, which 
would suppress all fine-structure lines due to an in¬ 
crease of dust competition for UV photons. This is 
further supported by the fact that the total gas cool¬ 
ing ([CII]+[OI]+CO+[CI]) decreases as a function of 
Lfir, meaning that dust becomes an even more efficient 
coolant as the infrared luminosity increases. 


Diaz-Santos 
star forma- 


The lack of a strong CO-line deficit shows that the 
bulk of the molecular gas heating is not affected by the 
mechanism suppressing the fine structure lines. This 
is interesting, since in a UV-photon heated environ¬ 
ment, suppressing the UV field implies a reduced heat¬ 
ing rate, and therefore also a lack of warm molecular 
gas might be expected. We also note that the integrated 
CO luminosity, that represents 10s of percent of the to¬ 
tal gas cooling, is much more than what is predicted by 
any pure PDR model, which give a CO cooling frac¬ 
tion of at most a few percent (3-5%) of the total gas 
cooling (e.g., Meijerink et al.|2011[ >. This result sug¬ 
gests that the CO luminosity may be powered by a dif¬ 
ferent heating mechanism, which does not lead to dis¬ 
sociation or ionization. 


4.3. Heating Mechanisms 

The CO molecule can be heated indirectly through 
(a) the photoelectric effect by ultraviolet (UV) pho¬ 
tons, (b) by fast electrons from directly ionized H and 
FF by X-rays, cosmic rays ( [Meijerink & Spaans 2005| , 
or (c) mechanical processes, which includes shocks 
and turbulence. X-rays heat gas by ionizing H and 
Hi directly, and these fast electrons then thermalize 
the molecular gas with an efficiency of 10%. UV pho¬ 


tons ionize PAHs and dust grains, and the resulting free 
electrons heat the molecular gas with a net efficiency of 
1-3%. In addition, the chemistry in an XDR is driven 
by X-ray photons instead of FUV photons that are able 
to penetrate deeper into the cloud without efficiently 
heating the dust at the same time. These X-rays are 
mostly produced by active galactic nuclei (AGN) or in 
areas of extreme massive star formation. Cosmic rays 
can also heat the gas by penetrating into cloud cen¬ 
ters, similarly to X-rays, and are typically produced 
by supernovae. Mechanical heating is another efficient 
source of gas heating. This is commonly attributed to 
turbulence in the ISM, which may be driven by su¬ 
pernovae, strong stellar winds, jets, galaxy mergers, 
cloud-cloud shocks, shear in the gaseous disk, or out¬ 
flows. 


To investigate the main mechanism heating the 
molecular gas, we use another diagnostic molecule, 
namely polycyclic aromatic hydrocarbons (PAHs). 
PAHs are carbonaceous, nanometer sized macro¬ 
molecules that contain 50-100 carbon atoms with an 
abundance of 10 7 per hydrogen atom (Tielens 2008 


The absorption of one far-UV photon is enough to heat 
the PAH molecule to a high temperature and will cause 
this molecule to emit in the characteristic bands at 3.3, 
6.2, 7.7, 8.6, and 11.2 /mi ( |Tielens|[2008| and refer¬ 
ences therein). Because PAHs are only fluorescently 
excited, and are easily destroyed by more energetic ra¬ 
diation, they are ideal tracers of UV heating. Thus, we 
can use the equivalent width of the 6.2 gm feature from 
the Spitzer Space Telescope (Stierwalt et al. 2013| , as a 
proxy for the UV energy density. In Figure [8j we com¬ 
pare the PAH equivalent width to the percentage of the 
total gas cooling done by CO, as calculated in Sec¬ 
tion 4.2 Class I objects have the steepest decreasing 


CO SLED, high PAH equivalent widths (EW), and low 
percentages of CO cooling. On the other hand. Class 
II and Class III objects with high percentages of CO 
cooling have low PAH equivalent widths. We note that 
the two objects with very low PAH equivalent widths 
are the most AGN-dominated objects in our sample 
(Mrk 231 and IRAS F05189-2524), where PAH de¬ 
struction by X-rays and nuclear hot dust combine to 
strongly lower the PAH equivalent width. The high 
CO cooling fractions of these objects can likewise be 
attributed to energy input by the AGN (through X-ray 
heating or mechanical energy input from an AGN- 
driven outflow). There are two galaxies that are Class 
I objects with lower PAH equivalent widths, and those 
are IC 1623 and NGC 7469. NGC 7469 is a starburst- 
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Fig. 8.— The PAH 6.2 /urn equivalent width from the 
Spitzer IRS (Stierwalt et al. 2013) > as a function of per¬ 
centage of CO cooling (for J > 4), as calculated in 
Section 4.2 Each galaxy is color-coded by class, red 
is Class I, green is Class II, and blue is Class III. In 
addition, we exclude Arp 220, NGC 4418, and Zw 
049.057 since their [OI] profiles are fully in absorption 
or show a complex line profile, making our CO cool¬ 
ing percentage inaccurate. Finally, we exclude NGC 
1365 since it is extended and we do not capture all of 
the CO emission in one SPIRE pointing, affecting the 
accuracy of our CO cooling percentage. 


ing galaxy with a Seyfert 1 nucleus and IC 1623 is a 
late-stage merger with a starburst nucleus. NGC 7469 
has a central starburst ring along with its AGN nucleus 
(e.g. Davies et al. 20071. This suggests that the Class 
I type CO emission is coming from the starburst ring 
that is encircling, but not directly affected by the AGN. 
These trends reinforce the idea that objects with a high 
CO cooling fraction, the CO is efficiently excited by 
something besides UV photons. 

We now compare our a parameter with the CO J=l- 
0 linewidth for the sources in our sample (Figure [9|. 
We correct the linewidth for inclination such that all 
galaxies are effectively turned edge on using the K- 
band axis ratio from 2MASS, with the exception of 
Arp 220 and NGC 6240. For Arp 220 we used the 
inclination derived by Scoville et al. (1997 1 using arc- 
second imaging of CO J=l-0 and for NGC 6240 we 
used the inclination derived by Engel et al. (2010) us¬ 
ing subarcsecond near-infrared imaging. The range of 
Class I objects is highlighted with a red. Class II with 
a green, and Class III with a blue background. 


Merging and interaction, molecular outflows, and 
random motions may all play a role in increasing 
the linewidths, but are not expected to dominate the 
linewidths of our target sources. The linewidth of 
CO J=l-0 is dominated by rotation, which is deter¬ 
mined by the mass of the central regions of the galaxy. 
Therefore, galaxies with high excitation (a) and high 
linewidths are more massive. However, the speed of 
rotation is also a proxy for the mechanical energy 
reservoir available in the galaxy nucleus since rotation 
leads to processes such as sheering and turbulence. 
If the fraction of mechanical energy that is converted 
into heating the molecular gas is constant, then high 
linewidth galaxies in Figure [9] would have more me¬ 
chanical energy available, and would be responsible 
for heating more of the molecular gas. However, es¬ 
tablishing this result would require a detailed study of 
the velocity fields of the molecular gas in all of our tar¬ 
gets, now possible with ALMA. Since such data is not 
available, thus we suggest that the linewidth provides 
an estimate of the available reservoir of mechanical en¬ 
ergy, but we remain agnostic as to the extent that this 
reservoir is actually tapped. 

Inspecting Figure [9] with these considerations in 
mind, the fact that low-excitation galaxies all have 
the smallest linewidths, while high-excitation galaxies 
have a range of linewidths has an attractive physical 
interpretation. In high excitation galaxies with high 
linewidths, mechanical energy input is a viable source 
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Fig. 9.— CO J=l-0 linewidth (FWHM) in km s“' 
plotted against a (Eq. [TJ. The different classes are 
highlighted with different background colors, red is 
Class I, green is Class II, blue is Class III. Most sam¬ 
ple galaxies are plotted with black dots, but some 
of those specifically mentioned in the discussion are 
highlighted. This also includes the three targets ad¬ 
dressed in Figures 1-3. 



Fig. 10.— The AGN contribution of the bolometric 
luminosity compared to the CO J=l-0 inclination cor¬ 
rected linewidth for Class III galaxies along with the 
Class II template galaxies presented in Figure [9] We 
have also included the rest of the Class II galaxies as 
transparent green points for reference. 


of excitation, while in high excitation galaxies with 
low linewidths, radiative energy input may be more 
important. Indeed, in Fig. 10, the lowest linewidths 
among high excitation galaxies are found in Mrk 231 
and NGC 4418 (an exposed and an obscured AGN, re¬ 
spectively). 


Class I and II objects with low linewidths are likely 
dominated by UV heating, since their CO ladders turn 
around somewhere before .1=7-6. For comparison, 
we have highlighted the position of our three exam¬ 
ple galaxies, NGC 7552 (Class I), Mrk 331 (Class II), 
and IRAS F17207-0014 (Class III), as well as some of 
the more famous galaxies in our sample. The major 
merger galaxies NGC 6240 and Arp 220 have a rela¬ 
tively high a values and linewidths around 650 krn/s. 
On the other hand, the strongest AGN in our sample, 
Mrk 231, lies in the low linewidth region of Class III, 
suggesting that its gas is radiatively heated, likely by 
X-rays from the AGN. This confirms the results of van 
der Werf et ah ( 201 Ob ), where they find the high-J 


CO excitation consistent with an XDR. The other two 
galaxies with AGN contribution are NGC 1068 and 
IRAS F05189-2524, both of which are Class II ob- 
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jects that lie in the low linewidth region, suggesting 
that they are also radiatively excited. However, since 
both of these objects lie in Class II (albeit on the border 
between Class II and III), it is unclear whether their ex¬ 
citation is from UV, higher energy photons, or cosmic 
rays. Both NGC 1068 and IRAS F05189-2524 have 
also been studied in detail by Spinogli o et al.| ( |2012| ) 
and Pereira-Santaella et al. ( |2014| > respectively. For 
NGC 1068, Spinoglio et al. ( 201 2j i and Garcia-Burillo 


|et al.| ( |2014 1 find that indeed, the excitation is due to 
either XDRs (in the circumnuclear disk) or PDRs (in 
the star-forming ring), which agree with our results. 
In IRAS F05189-2524, [Pereira-Santaella et al.ld20T4| 
find that there is a large contribution from mechanical 
heating in this source, as traced by a shallow H 2 tem¬ 
perature distribution, yet they cannot rule out the AGN 
as an important heating source for the molecular gas. 

In order to examine this method vis-a-vis the role 
of AGNs and to determine additional heating sources, 
and to check if there are any underlying biases in the 
J=l-0 linewidth, we can compare the inclination cor¬ 
rected linewidths with the percentage AGN contribu¬ 
tion for each galaxy. In the case of very disturbed ma¬ 
jor mergers, such as Arp 220, Arp 299, IRAS 17208- 
0014, and NGC 6240, the inclination is difficult to 
determine, and thus the corrected linewidth can be 
over- or underestimated. The uncorrected linewidths 
are shown in Table 1 for comparison. The AGN con¬ 
tribution can be estimated using both the 15/30 yum 
flux density ratio (/15//30)and the [NeV]/[NeII] ra¬ 


tio along with the prescriptions from Veilleux et al. 


(2009). We take the /15//30 ratio from Stierwalt et al. 
(2013|> and the [NeV]/[NeII] ratio from Inami et al. 


H3). We average the results of the AGN contribu¬ 
tion from both methods and find average AGN contri¬ 
butions ranging from to 0%-95% of the bolometric lu¬ 
minosity. For the high excitation sources (Class III) we 
use the average AGN contribution in combination with 
the inclination corrected linewidth to separate mechan¬ 
ical heating from AGN heating. In Figure [TO] we plot 
the AGN contribution on the x-axis and the inclination 
corrected CO J=1-0 linewidth along the y-axis. This 
scatter plot shows that for our high-excitation galaxies, 
very high linewidths are not associated with high AGN 
contributions and conversely the galaxies with high 
AGN contributions do not display high linewidths. We 
also plot the Class II galaxies as green points (the tem¬ 
plate galaxies are marked) for comparison. Most Class 
II galaxies have low AGN contributions, while they 
have a large range of linewidths. 


For example, we can compare Arp 220 and Mrk 
231. Both galaxies have high a values, but Arp 220 has 
a high linewidth and a low AGN contribution, while 
Mrk 231 has a low linewidth and a high AGN contri¬ 
bution. Comparing NGC 6240 and Arp 220 reveals 
that although both have high linewidths and a values, 
NGC 6240 also has a high AGN contribution. This 
suggests that both the AGN and mechanical processes 
are contributing to heating the gas in NGC 6240. 


It is difficult to conclude anything definitive about 
objects that show average or typical values of a or 
linewidths. ft is also important to note that although 
there is a non-negligible contribution of heating from 
mechanical processes or the AGN in Class Iff galax¬ 
ies, the gas is still mostly heated through UV heating 
processes. We caution the use of any two of these di¬ 
agnostics alone, since for example, a Class Iff object 
with a low linewidth may still be mechanically heated 
by small-scale turbulence that would not produce an 
observable global line broadening effect. In this case, 
the AGN contribution would be low, but the a would 
be high, discounting AGN heating. Using all three 
parameters simultaneously allows for a qualitative es¬ 
timate of which additional processes are exciting the 
warmest molecular gas. In order to fully understand 
the excitation mechanisms and physical parameters of 
the molecular gas, an additional detailed modeling of 
the 12 CO, 13 CO, and dense gas tracers (HCN, HNC, 
HCO + , etc.) is required (e.g. Rosenberg et al. 2014a 1 . 


5. Conclusion 


We report the initial results of the Herschel Open 
Time Key Project Herschel Comprehensive (U)LIRG 
Survey (HerCULES). Both Herschel/SPIRE spectra 
and Herschel/PACS [OI]63, [OI] 145 and [CII] line 
profiles of a sample of 29 galaxies spanning an order of 
magnitude of infrared luminosity were analyzed. Our 
main results are summarized below: 


• We separate our sample of luminous galaxies 
into three qualitative classes based on the shape 
of their CO ladder, quantized by the parameter 
a, which we define as the ratio of the high-J to 
mid-J CO transitions. 

a = (LcO j =11-H)+LcO j =12-U+LcOj=13-12)/(LcOj=5-4 + 
LcOj= 6-5 + L C Oj=l-(,) 

Class I (a < 0.33) is characterized with a CO 
SLED peak around J =5-4 and a steep decline to¬ 
wards higher J transitions. Class II (0.33 < a < 
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0.66) has a CO SLED peak around .1=7-6 and 
a shallower decline towards higher J transitions. 
Class III (a > 0.66) shows very flat CO ladders. 
We present the spectra of three exemplary galax¬ 
ies for each of these three categories. 

• We find that molecular gas excitation (approxi¬ 
mated by a) is well correlated with the infrared 
color (as traced by S 60 /S 100 ). and not as well 
correlated with the Lir. 

• The cooling budgets of the galaxies are pre¬ 
sented. We find that the percentage of cool¬ 
ing from each species ([CII], [Sill], [OI], and 
[Cl]) appears to be constant over the full range 
of Lfir. There is indication of a slight increase 
in the percentage of CO cooling at higher Lfir. 

• We find [CII]/FIR, [SiII]/FIR, and [CI]/FIR ra¬ 
tios consistent with deficits (i.e., the [CII]/FIR 
ratio decreases with increasing L fir), and also a 
weak [OI] deficit for the high luminosity galax¬ 
ies. On the other hand, we observe no CO 
deficit, the CO/FIR ratio is very constant for all 
L F ir, with the exception of NGC 6240. Thus, 
the fine structure line deficits reflect a decrease 
in the total gas heating efficiency with increasing 
L fiR. The fact that we observe a deficit in all fine 
structure lines but not in the molecular gas sug¬ 
gests that the mechanism responsible for heating 
the [CII], [OI], [Sill], and [Cl], is not the same 
mechanism responsible for heating the CO. The 
CO may instead be affected by a heating mech¬ 
anism that is immune to this deficit. We also 
find that the total neutral gas cooling per FIR 
decreases as a function of h FIR . 

• Using the PAH 6.2 /mi equivalent width as a 
proxy for the importance of massive star forma¬ 
tion and therefore UV excitation, we find that 
when the cooling efficiency of CO is high, the 
amount of UV heating is low. This again in¬ 
dicates that CO is more efficiently heated by a 
mechanism not directly related to UV radiation. 

• We suggest a qualitative schematic based on a, 
the CO J=l-0 linewidth, and the AGN contribu¬ 
tion, that helps indicate which additional mech¬ 
anism, if any, is heating the gas. Class I galax¬ 
ies with a low a (a < 0.33) do not require 
any heating in addition to UV-heating to explain 
the observations. Class III galaxies with high 


linewidths and low AGN contributions probably 
require mechanical, in addition to UV heating. 
Class III galaxies with narrow linewidths and 
large AGN contributions are experiencing exci¬ 
tation from harder radiation (X-rays or cosmic 
rays). Class III objects with wide linewidths and 
high AGN contributions are composite galaxies 
that are being heated by both mechanical pro¬ 
cesses and the AGN. For the objects that have 
median a, linewidth, or AGN contribution val¬ 
ues, such as many Class II objects, it is not pos¬ 
sible to discriminate which heating mechanisms 
are affecting the gas without additional informa¬ 
tion. 
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